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Abstract 

The solitons of affine Toda field theory are related to the spin- 
generalised Ruijsenaars- Schneider (or relativistic Calogero-Moser) mod- 
els. This provides the sought after extension of the correspondence be- 
tween the sine-Gordon solitons and the Ruijsenaars- Schneider model. 



1 Introduction 



The purpose of this letter is to relate the solitons of the affine Toda system to 
the spin extensions of the Ruijsenaars- Schneider model, the latter also being 
known as the relativistic Calogero-Moser model. This work generalises the 
known connection between the solitons of the sine-Gordon equation and the 
(non-spin) Ruijsenaars-Schneider model. The connection made here should 
be viewed as part of a larger programme that seeks to identify classes of solu- 
tions of PDEs with finite dimensional mechanical systems, whereby the evo- 
lution of the solutions to the PDE is expressed as a dynamical system on the 
(finite-dimensional) moduli space of solutions. Thus, for example, the pole 
solutions of the KP equation |], || and its reductions (such as KdV 0) are 
related to the non-relativistic Calogero-Moser model, while the sine-Gordon 
solitons are related to its relativistic counterpart, the Ruijsenaars-Schneider 
model (see f|, ||). This programme also extends to include the peakon solu- 
tions appearing in fibre-optics and shallow water waves which have associated 
mechanical systems @, 0|. A similar connection may well underlie the ap- 
pearance of finite dimensional mechanical systems in the study of various 
models possessing duality @, ||. As the dynamics of mechanical systems are 
often easier to understand (or simulate) than the equations of motion for a 
field theory, such a programme aims at giving qualitative information about 
field theories by an appropriate reduction of degrees of freedom. The recent 
work of Babelon, Bernard and Smirnov [|l(J may be viewed as taking this 
correspondence between field theories and mechanical systems beyond the 
classical to include the quantum regime as well, though the ability to focus 
attention solely on a fixed N-particle sector of the full quantum Hilbert space 
appears to depend crucially on the model. Our work will reveal further new 
features in such correspondences, as well as provide a sought after general- 
isation of known results about the sine-Gordon model to the case of affine 
Toda solitons. 

Ruijsenaars and Schneider's seminal work || showed that the soliton 
solutions of a variety of equations were related to dynamics built from the 
Hamiltonians (with canonically conjugate variables qj, pj) 



N 



N 




J 



and equations of motion (for either H±) 




(1) 
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An appropriate scaling limit of H = H + + H_ yields a system of Calogero- 
Moser type. In particular, the eigenvalues % e qj of an N x N matrix associated 
with the tau function describing an N-soliton solution of the sine-Gordon 
equation evolve according to flU). The q/s and p/s may, at least when they 
are well separated, be related to the positions and rapidities of N constituent 
single-solitons; the dynamics of the system encodes the various soliton phase 
shifts. (More details of this will be given below.) Thus the system governed 
by H + describes how the space-time trajectories of the 'constituent' solitons 
interact. Of course the same system may be described via the inverse scat- 
tering transform by a free system with linearly evolving data: the point of 
the Ruijsenaars-Schneider description is to make greater contact with the 
particle description of the soliton. 

Viewing the sine-Gordon model as the afhne Toda system (with imag- 
inary coupling) a natural question to ask is how the above results generalise 
to other affine Toda systems. These systems have been extensively studied 
in recent years both for real and imaginary couplings. In the real coupling 
regime a beautiful structure was uncovered and exact S-matrices have been 



conjectured for the theories (see [IT[ [12], [13| and references therein). The 



imaginary coupling regime has also been investigated and classically the soli- 



tons have real energy- momentum although the Lagrangian is complex ([[14 
and references therein). Spence and Underwood [15 have recently used this 
work to obtain the symplectic form on the space of affine Toda solitons but 
a dynamical description generalising the sine-Gordon/Ruijsenaars-Schneider 
correspondence has proved elusive. The purpose of the present letter is to 
give this generalisation. Just as the affine Toda systems generalise the sine- 
Gordon model, there are spin-generalisations of the Ruijsenaars-Schneider 
systems, and it is these systems which describe the dynamics of the affine 
Toda solitons. These models (which have been most studied in the A n set- 
ting) are the relativistic extension of Gibbon and Hermsen's spin generali- 
sation of the original Calogero-Moser model [19. One new feature we have 



found in our correspondence is the appearance of new degrees of freedom, 
the internal spins of the model. Although not needed to describe the solitons 
of the affine Toda system, these spins determine the matrix that diagonalises 
the Lax pair. We will comment further on this later in the letter. 

An outline of the letter is as follows. First we will review the construc- 
tion of affine Toda solitons, and then in section 3 consider the reduced sym- 
plectic form of the theory. We are then in a position to relate the affine 
Toda solitons to the spin-generalised Ruijsenaars-Schneider model in section 
4. For the purposes of this letter we shall limit our discussion to the An 
case, both for simplicity and to make clear the generalisation of the sine- 
Gordon/ Ruij senaar s- S chneider corresp ondence . 
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2 The A$ Affine Toda Solitons 

For the An affine Toda theory with imaginary coupling, the equations of 
motion read 

2 

d+d-fa + ^_( e ^fe-^+x) _ e ifWj-i-*i)) = o, (2) 

j = 0, l,...,n. Here ± denotes differentiation with respect to light-cone 
coordinates x± = -^(t ± x), and the indices on the components of the field 
<p are read modulo (n + 1) where necessary. We shall be considering the 
solitonic sector of the theory, which means assuming ^™ =0 <f>j = ( m other 
words, discarding the free field part of 0). 

There are various ways to construct and parametrise soliton solutions to 
(fj). Perhaps the simplest methods to implement from a practical point of 
view are the application of the Backlund transformation derived by Fordy 
and Gibbons jI7| or the bilinear formalism developed by Hirota [pl|. There 



are also the powerful vertex operator techniques which make full use of the 
representation theory of the An algebra [[14]] . While the latter approach is 
currently the most popular, we wish to make contact with the original work 
of Ruijsenaars and Schneider [0], which made much reference to the soliton 
formulae of Hirota. Hence we choose to start from the form of the N-soliton 
solution of (H) derived by Hollowood via Hirota's direct method. The 
l-th component of the field is given by 

e iM = Zhl (3) 

n 

where the tau function n is of the form 

T l = Yl eXP ( e 3 e kB jk + e 3Q,l( X +> X -) ) ■ ( 4 ) 
e \j<k j J 

In the above the e indicates a summation over all possible combinations of 
6j taking the values or 1, and the indices j and k take values in {1, iV}. 
We will explain shortly what the various terms in (f|) mean, and how we have 
parametrised the An affine Toda solitons. For the moment we would like 
to comment that expression (Q) is a rather generic form of the soliton tau 
function for an integrable PDE, the precise nature of Bjk and depending 
on the particular PDE being considered; it may be viewed as a degeneration 
of the theta function solutions of the PDE given via algebraic geometry 
in which the e/s run over all of the integers. Ruijsenaars and Schneider 
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succeeded in making the connection between their relativistic Calogero-Moser 
systems and soliton solutions of the sine-Gordon and KdV equations, among 
others, by showing a direct correspondence between the coordinates of the N- 
particle system and the parameters of the N-soliton solution. An important 
part of the correspondence was that all the tau functions of form (f|) being 
considered in || could be written in terms of determinants like 



for suitable matrices M. In what follows we express all the iV-soliton solu- 
tions of the An^ affine Toda theory in this way, and thereby obtain a relation 
to spin-generalised Ruijsenaars-Schneider systems. 

First of all we should explain the parameters of the Toda N-soliton which 
appear in (§]). Each soliton has a rapidity denoted by rjj, a position param- 
eter denoted by aj, and a discrete parameter 9j taking values in {|rf|A; = 
1,2, ...,n} (so that exp(i6j) is an (n + l)th root of unity). The rapidities 
are all real, while the aj are pure imaginary for solitons (there are different 
reality conditions for other types of solution e.g. breathers). The different 
values of 6j give n different species of soliton in the An affine Toda theory 
whose masses are 2msin(^/2). We also need to define 



With this choice of parameters, the terms in the sum (|]) are given by 



(To make a comparison with the vertex operator formulae, we note that in 
terms of the notation of reference [[OJ, we have Bjk = log(Xj^), aj = Qj. 
We will deal with the general formalism elsewhere.) 

We are now ready to write the tau functions as determinants. In fact 
Olive, Turok and Liao |20] found that determinants naturally arose when 
they derived the N-soliton solution by the Backlund transformation, but the 
matrices involved are not of the right form for our purposes. Instead we set 
Xj = aj(fj,'j — fij) exp (v / 2m(e~' ,:, i + — e Vi xJ) sin(0j/2)), and define N-by-N 
matrices V , 6 by 



det (1 + M) 



l*f = expfo- ± -i9j). 






and 

e = diag(0 1 ,0 a ,...,0 JV ). 

Then we find that 

n = det (l + e iie/2 Ve m/2 ). (6) 

To verify (^) it is necessary to expand the determinant on the right-hand 
side in terms of the principal cofactors of V, and then use Cauchy's identity: 

dct V 1 \ _ TT 1 tt (^ + -^)(^ 7 -^) 

W - ^ / t $ ~ ^ i<k W ~ Vk)(vj - f4) ' 

Writing the right-hand side of (f|) in terms of these new parameters and 
comparing with the cofactor expansion gives the required result. Note that 
in the Af> (sine-Gordon) case the 9j must all take the value tt, which means 
that the matrix exponentials appearing in ffl) are multiples of the identity, 
and we reproduce the standard result 

m = p -m, = det ( x ~ v ) 

det (1 + V)' 

3 The Reduced Symplectic Form 

In this section we describe the phase space of the iV-soliton solution in terms 
of its symplectic form, before describing how spin-generalised Ruijsenaars- 
Schneider systems arise in the following section. The phase space of the affine 
Toda system has the standard symplectic form 

/oo 
(6(pt A, 5<p) dx. (7) 
-oo 

On substitution of the iV-soliton solution into (|7j), one obtains (after an 
integration) the reduced symplectic form on the iV-soliton phase space. In 
practice it is not possible to perform the integration for anything other than 
the one-soliton solution |15| (except for the sine-Gordon case, where Babelon 
and Bernard succeeded in showing that the integrand could be written as an 
exact derivative for both the one- and two-soliton [|j]). For the one-soliton 
phase space, the reduced symplectic form is (up to an irrelevant numerical 
factor independent of 9) 

w (i) = _^ a drj. 
a 

The intractability of the integral (^) for the general iV-soliton solution 
does not matter, as it is a standard result that as t — > ±oo (the out /in 
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limits) the iV-soliton decomposes into a superposition of iV one solitons with 
a shift of the parameters. So the symplectic form may just be written 

- <K) = £^Arf,r = £^A*?. (8) 

j i j 3 

By direct calculation using the formula (0) for the tau functions of the N- 
soliton solution, we find the relations between the out /in parameters and the 
standard ones: 

r.in _ out _ 

(and similarly for a° ut with the inequality reversed). This agrees with the 



formulae of Spence and Underwood |L5[ obtained via vertex operator argu- 
ments, where in their notation dj = Qj and exp(Bjk) = Xj jk . So substituting 
for the in parameters into (pj), we obtain the N-soliton symplectic form as 

U {N) = ^ -^1 a drjj + ^2 E ik{T]) sinhfo - rj^drjj A dn k , (9) 

where 



cosh(?7 j - r) k ) - cos((6' i - 6 k )/2) cosh^ - n k ) - cos((^ + k )/2) ' 



We observe that W M is clearly real if we choose the rjj to be real and the 
dj to be pure imaginary (which in the A± case coincides with the condition 
on dj for sine-Gordon solitons given in Q). This means that the matrix V 
defined in ([5]) is anti-hermitian, which will be important in the next section 
when we look at the dynamics of the eigenvalues of V. 



4 Ruijsenaars-Schneider Systems 

Here we consider how the eigenvalues of the matrix V evolve with respect to 
each of the light cone coordinates, and find that spin-generalised Ruijsenaars- 
Schneider equations result. Since V is anti-hermitian, it may be diagonalised 
with a unitary matrix U: 

Q := UVU ] = diag(iexp(gi), . . . , i exp(q N )) . 
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If we let a dot denote then V satisfies 

V=±{AV + VA), (10) 

for the constant diagonal matrix A given by 

A = diag(±\/2mexp(=F77i) sin(6>i/2), ±v / 2mexp(=F?77v) sm(6 N /2)). 

Now let Uj denote the jth row of U (considered as a column vector, so that 
the Uj are the left eigenvectors of V). Define the Lax matrix L by 

L jk = u[Auj, 

so that 

L = UAU j . 
Then L satisfies the Lax equation 

L = [M,L], (11) 

for M = UW . Differentiating the definition of Q gives 

Q = [M,Q] + UVrf, 

and after substituting for V from (|10| ) and using the definition of L we find 
the identity 

LQ + QL = 2{Q+[Q,M\). (12) 
Upon setting Qj = i e 9j , then ([l!|) in components reads 

Ljk(Qj + Qk) = 2(Qj5jfc + Mjk(Qj — Qk)), 

which yields 

Ljj = qj (13) 

and 



1 ( Qj + Qk 



) L 3k = \ coth((g i - q k )/2)Lj 



(for j ^ k). When V is diagonalised we may always choose the phases of the 
left eigenvectors Uj so that Mjj = Ujiij = 0. Substituting these into the Lax 
equation (jTTf) produces the equations of motion: 

Aw = Qj = ^ coth((g j - q k )/2)L jk L kj , (14) 



Ljk = \ coth((# - q k )/2)(q k - qj)L jk , , 

+ J2 ¥j , k |(coth((g, - qi )/2) - coth((g ; - q k ) /2))L 3l L lk 1 ) 

(j 7^ fc). These are in fact the spin-generalised Ruijsenaars- Schneider equa- 
tions with certain constraints, although to see this requires comparison with 



the formulae of Krichever and Zabrodin 21 



In [^TJ the generalised Ruijsenaars- Schneider model is defined in terms of 



N particle positions Xj and their internal degrees of freedom (spins) given 
by /-dimensional vectors a 3 - and /-dimensional covectors fej, subject to the 
equations of motion 

Xj = ^{b)a k )(b\a,j){V{xj - x k ) - V(x k - Xj)), (16) 

k+j 



y]a h (6^0j)V(gj -g fc ), (17) 

k+j 

b\ = -Y, h W^k)V{x k -Xj). (18) 

k+j 

The potential V is expressed in terms of the Weierstrass zeta function or its 
rational or hyperbolic limits. To make contact with our equations we set 
Xj = qj and choose the hyperbolic potential 



V(<?j -Qk) = \ coth((^- - q k )/2). 



Then ( |T6| ) generalises ([!]). In [^TJ the spin degrees of freedom were real, but 
here we allow them to be complex, and identify them with the eigenvectors 
of V by setting 

fet =Up dj = Auj. 

So we have taken / = N, and in fact our spins are expressed entirely in terms 
of the eigenvectors of V and the constant matrix A; in particular the fot must 



form an orthonormal basis. In the notation of the components of the 
Lax matrix are given by 



L jk - b k 



k a j- 



There are various other constraints that we have imposed on our system. 
First the equations (|16|), ([17]) and fll8l) have the scaling symmetry 



a j -> «i a i> h ) — &}- 

(Xj 



The corresponding integrals of motion are Xj — bjCij, and setting them to zero 
and rewriting them in terms of our coordinates shows that this is equivalent to 
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equation ([13]). Similarly our requirement that Mjj = is another constraint 
on the system. Now given these constraints we find that from the definition 
of L in terms of the spins we can compute L^. So for j = k (jTj|) is equivalent 
to (HD, while for j ^ k ($7$) and (0) yield (0). 

To make the correspondence between the solitons and the many-body 
system clearer, it is worth considering the sine-Gordon case in more detail 
and comparing it with the general situation. The results about sine-Gordon 
solitons are explained in detail in [|J , and we have kept our notation as similar 
to this reference as possible to make comparison easier. The first thing to 
observe is that in the case only knowledge of the qj is required to specify 
the field components, as we have 



N 

e i/30o = e -Wi = 

3=1 



1 — i exp(qj 
1 + i exp(qj 



In the general case the presence of the matrix e de>//2 in the expression for the 
tau functions (g|) means that knowledge of both the spin vectors Uj (which 
make up the matrix U) and the qj is required to evaluate these determi- 
nants. The essential difference is that for sine-Gordon there is only one soli- 
ton species, while in the An case there are n different species corresponding 
to the different allowed values of 6j. This difference is also apparent at the 
level of the equations of motion. In fact when we differentiate the matrix V, 
in the case of sine-Gordon we find from (ITOl) that 



V = i(ee j ) 

for a certain vector e. But then conjugating the equation (|T0| ) with U we 
obtain 

ie$ = ^(LQ + QL), 

where e = Ue. Actually e is a real vector, and in terms of its components ej 
we have 

U 



_ 2 e j e k 

exp(r/ ; ) • e.xpO//,)' 



Since we know the diagonal elements of L explicitly in terms of the qj (from 
fljjD ) the above formula means that we then know all the e.j and hence the 
off-diagonal elements of L are found to be 



Ljk 



cosh((g i - q k )/2)' 

This may then be substituted into (|14]),(|T5|) to give the ordinary (non-spin) 
Ruijsenaars-Schneider equations. In this case flH[ ) yields ([I]) and (|TS|) is a 
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consequence. Babelon and Bernard have shown [Q that there is a canonical 
transformation between the soliton parameters and the dynamical variables 
qj , (jj (more precisely, they formulate this in terms of the variables Qj = 
i exp(qj)). We discuss how this could possibly be extended to the An case 
in our Conclusion. 

5 Conclusion 

We have shown the connection between spin-generalised Ruijsenaars-Schneider 
systems and affine Toda solitons. The soliton tau functions are deter- 
mined by the positions qj of n particles on the line as well as an orthonormal 
set of n-dimensional spin vectors Uj, which are together subject to the equa- 
tions of a constrained spin-generalised Ruijsenaars-Schneider model. This 
extends the known result for the sine-Gordon equation, where the spins are 
no longer part of the dynamics and there is a canonical transformation be- 
tween the positions and momenta of the particles and the parameters of the 
solitons. For the general case such a transformation is no longer apparent, 
although we note that the iV-soliton phase space is still of dimension 2N, 
and so it is worth exploring exactly how the extra spin degrees of freedom 
are absorbed in the transition from the dynamical variables to the soliton 
parameters. Also it would be interesting to see what role the spins might 
play in the quantum theory. Finally there remains the extension to the other 
affine algebras and elucidating the connections to the vertex operator con- 
structions mentioned at various points in the text. We intend to pursue these 
points in the future. 
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